o 

(N 
O 

Q 



C/2 ■ 

ctf ■ 

i ' 
^— > ■ 

C ■ 

c3 ■ 

=5 

■ 

e 

C 

o 
o 



> 

in 

(N 

in 
in 



X 



Typeset with jpsj3.cls <ver.l.l> 

One- Dimensional Integrable Spinor 
BECs Mapped to Matrix Nonlinear 
Schrodinger Equation and Solution of 
Bogoliubov Equation in These Systems 

Daisuke A. Takahashi* 

Department of Basic Science, the University of Tokyo, 
Tokyo 153-8902 

KEYWORDS: matrix nonlinear Schrodinger equation, spinor 
Bose-Einstein condensate, Bogoliubov equa- 
tion, squared Jost function. 

Nonlinear Schrodinger equation(NLSE) is one of the 
oldest integrable nonlinear equations, which was solved 
by the inverse scattering method(ISM). 1 ) In the con- 
text of condensed matter physics, NLSE is also called 
Gross-Pitaevskii equation, and it describes the dynamics 
of one-dimensional(lD) Bose-Einstein condensate(BEC). 
Recently, integrable ID BECs with internal spin degree 
of freedom (spinor BECs) have been discovered. 2, 3 ) First, 
the spin-1 integrable system has been discovered 2 -* by 
finding a mapping to the matrix NLSE(MNLSE): 4 ) 

-d 2 x Q±2QQ^Q. 



id t Q 



(1) 



Here Q is a square matrix, and the minus(plus) sign cor- 
responds to the self- (de) focusing case. In the context of 
BEC, the minus(plus) sign represents the system of at- 
tractive(repulsive) bosons. Subsequently, integrable BEC 
systems are found for every integer spin-n. 3 ) The energy 
functional of integrable spinor BECs are given by 3 ) 



H 



dx 



E 



m— — n 



dvb Tl 



dx 



± 



lei 



(2) 



where ip m is a condensate wavefunction with magnetic 
quantum number to, p — J^ m __ n \ip m \ 2 is a number 
density, and 6 = EL-nt" 1 )™^™'/'-™ is a singlet pan- 
amplitude. The coefficient 2™ _1 can be eliminated by a 
choice of normalization, but we keep it for later conve- 
nience. 

In this short note, we show that the above integrable 
systems (2) are all mapped to MNLSE, and solve the 
Bogoliubov equation of these systems. While the integra- 
bility itself has been already proved by using a different 
Lax pair, 3 ) the mapping to MNLSE has the following 
advantages: (i) The explicit expression for 7V-soliton so- 
lution is already given under both vanishing 4 ) and non- 
vanishing 5 ) boundary conditions, (ii) With the aid of the 
theory of squared eigenfunction(or squared Jost func- 
tion), 6,7 ) solutions of Bogoliubov equation can be ob- 
tained. Since there exist various kinds of low-energy ex- 
citations in spinor BEC systems, 8 ' it would be an in- 
teresting future work to solve the scattering problem of 
low-energy excitations in the presence of an external po- 
tential 9 ) by using the solution given in this short note. 

Mapping to MNLSE — We first write the mapping ma- 
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trix and its (anti)symmetrization, and next explain how 
to find it. The mapping from spin-n integrable spinor 
BEC to 2™-dimensional MNLSE is given by 



Vn=l 



(3) 



Here <7j(i = x,y,z) are Pauli matrices, o~± = ~^(?x ± 
icr y ), I is 2 x 2 identity matrix, and A® 1 := A <g> • • • (g> A. 



i times 



(for i — 0, we define ^4®° := 1.) The energy functional 
which yields Eq. (1) is 

/0Q f dQ 



H 



d.r 



tr 



\ dx dx 



±tr(QtQQtg) 



(4) 



Substituting the matrix (3) to this functional, one can 
obtain the functional (2). 

Furthermore, we can (anti)symmetrizc the above ma- 
trix by the unitary matrix 



V = 



(n : even) 



(ia y <g> cr^)®^" 1 )/ 2 ® (ia y ) (n : odd) 



(5) 



The matrix Q = QV then satisfies Q T = (-1)5™("- 1 )Q. 

We briefly summarize how to construct the matrix (3). 
In the theory of Lie algebra, 10 ) writing the elements of 
Cartan subalgebra as Hi , . . . , H r (r: rank of this algebra) 
and raising or lowering operators as E a (which changes 
weight from p to p + a), an irreducible representation D 
is determined from the following relations: 

Hi D) = p t D) , (6) 

E a \p,D)=N a ^\p + a,D). (7) 

Here p = (pi, . . . ,p r ) is a weight vector, and N at ^ is a 
representation-dependent constant. The irreducible ten- 
sor operator 11 ) corresponding to this representation, D, 
is characterized by the following commutation relations: 



[Hi, T^d] = Pi T^.d , 



(8) 
(9) 



Now, let us consider the algebra so(2n+l). We con- 
struct the tensor operators corresponding to 2n + 1- 
dimensional fundamental representation (we write it as 
2n + l — it is unitary equivalent to the definition of the 
algebra itself.). As for the matrix representation of gen- 
erators, His and E a s, we use the 2™-dimensional spinor 
representation. Therefore tensor operators also become 
2 n -dimensional matrices. In this representation, genera- 
tors are represented by the n-fold tensor product of Pauli 
matrices. 10 ) One can determine the matrix elements of 
T^,2n+i (m = 0, ±ei, . . . , ±e„) through the commutation 
relations (8) and (9), and obtain the matrices appearing 
in Eq. (3). 

Bogoliubov equation — Next, we solve the Bogoliubov 
equation of these integrable systems. Since we are in- 
terested in an application to the scattering problem of 
collective excitations, 9 ) we discuss the problem of non- 
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vanishing boundary condition, 5 ' and consider MNLSE of 
the sclf-defocusing case with a chemical potential term: 



id t Q = -{d 2 x + n)Q + 2QQ^Q. 



(10) 



The Bogoliubov equation, though its standard deriva- 
tion is diagonalization of second-quantized Hamiltonian 
in mean field approximation, is easily obtained by substi- 
tuting Q = Q + SQ in Eq. (10) and discarding the higher 
order terms of 5Q. Rewriting (SQ, — SQ^) — (U,V), one 
obtains 

id t U = -(dl + n)U + 2(QQ^U+UQ f Q-QVQ), 
id t V=(d 2 x +n)V-2(Q^QV + VQQ^-Q f UQl). ^ 

ISM of MNLSE is formulated through the following ex- 
tended Zakharov-Shabat type eigenvalue problem: 4,5 ' 

JL(f\ = (-& Q\(f 

dx\g) W iVVfl 



(12) 



Iff 
dt\g 



2AQt _ i^Qt 



i(-2A 



2AQ + \d x Q \(f 
i(2X 2 -^+Q^Q)J\g 

(13) 



Here A is a spectral parameter, and / and <?, which 
are called Jost function, are n-component vectors when 
Q is an n x n matrix. The compatibility condition 
d x d t (f,g) T = d t d x (f,g) T reproduces Eq. (10). 

Assume that Q is (anti)symmetric: Q = eQ T (e = ±1). 
We can then show by a straightforward calculation that 
if (/i,.gi) T and (f 2 ,g 2 ) T are the solutions of Eqs. (12) 
and (13) with the same A, (U, V) = (fif£, —egig^) sat- 
isfies Eqs. (11). Thus, the squared Jost function gives a 
solution of Bogoliubov equation. 

Example — As an example, let us consider the spin-1 
case in the presence of one dark soliton. 12 ' We want the 
stationary solution of Bogoliubov equation with eigenen- 
ergy e, which is obtained by the replacement id t — > e in 
Eqs. (11). Correspondingly, we must consider Eqs. (12) 
and (13) with idt — > e/2, since a quadratic form of / and 
g gives a solution of Bogoliubov equation. 

In order to simplify the mathematical descriptions, we 
first prepare several notations for the solutions of one- 
component NLSE (i.e., a scalar BEC). The one dark soli- 
ton solution in the comoving frame is 

(f) s {x;x ) = -e^e ipx ]p + iqt<mh{q{x-x ))\ , (14) 

where p = — A cos ip, q = A sin ip, A > and ip e (0, 7r). 
The subscript "s" means scalar. For this solution, the 
chemical potential becomes p, — p 2 + 2Ag. We also de- 
fine 4> s (x;— oo) := Xoe lpx . The solution of Bogoliubov 
equation is given by (u, v) = (J 2 , — g 2 ). Considering the 
solutions in the uniform system, one can show that the 
wavenumber of an excitation k and the spectral parame- 
ter A are related as k = — ej (2 A — p). The expression for 
(/, g) in the presence of one dark soliton is given by 7 ' 



f s {x;x ) 
gs(x;x ) 



e ikx/2 x 



(15) 



' e^P x +^/ 2 [iqt&nh(q(x-x )) + (k/2) + (e/2k)] , 
^-i e -Kpx+v)fc[i q t a nh(q(x-x )) + (jfe/2) - {e/2k)])' 



where the wavenumber k satisfies the dispersion re- 
lation (e - 2kp) 2 = k 2 {k 2 + 4A§). We also define 
/ B (a;;-oo) := c^+^+^Z 2 [iq + (A/2) + {e/2k)\ and 
g s {x; -oo) := ^ c A{k- P )x- v \/2 ^ + _ ( e / 2 fc)]. 

Having prepared the notations, let us move on to 
the spin-1 system. The symmetrized matrix is given 
by Q — ( , ^} r- ^°/^V and correspondingly, Bogoli- 

ubov wavefunctions become U = ( "/ u °/^ ) and 

V M o/V2 u_i / 

V = ( v ' Vo/V2 ~ ). The one dark soliton solution 12 ' can 

V"Uo/v2 v-i J 

be diagonalized as 



Q = D 



(j) s (x;xi) 






(f> s (x;x 2 . 



D 1 



(16) 



Here D is a rotation matrix of spin- 1/2. (Note that the 
rightmost term is not but D T . 12 ') If either X\ or x 2 
equals to — oo, Eq. (16) represents a ferromagnetic soli- 
ton. 12 ' If both arc finite, it represents a polar soliton. 12 ' 
Jost functions are given by 



(i = 1,2), 



(17) 



'A = f f s {x;Xi)Dei 
vS/ \g s {x\Xi)D*e, 

where e\ = (1,0) T and e 2 = (0, 1) T . From these Jost 
functions, we can construct the solutions of the Bogoli- 
ubov equation as 

'U\ f f s (x; Xl ) 2 D(I + a z )D T 
V) \-g a {x; Xl ) 2 D*{I + <j z )D\ 



(18) 



f a {x;x 2 ) 2 D{I-a z )D T 
-g a {x-x 2 ) 2 D*{I-a z )D^ 

f s (x; xi)f s (x; x 2 )Da x D T 
-g s (x;xi)g s (x;x 2 )D*cr x D J< 

In the integrable case, the dispersion relations of excita- 
tions of the spin fluctuation and the density fluctuation 
are degenerate. (See Ref. 8 for the general case; Note that 
the asymptotic form of Eq. (16) is always polar state, re- 
gardless of whether the soliton is ferromagnetic or polar.) 
Therefore any linear combination of the above Eq. (18) 
becomes the solution. However, if one wants to separate 
the spin fluctuation part and the density fluctuation part, 
one needs to make a special linear combination of them. 

Note added: After submission, we became aware of 
Ref. 13, in which the reductions of MNLSE equivalent 
to ours are derived. The author thanks T. Tsuchida for 
informing him of this work. 
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